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A system of three spin-1/2 atoms allows the construction of a reference-frame- free (RFF) qubit 
in the subspace with total angular momentum j = 1/2. The RFF qubit stays coherent perfectly as 
long as the spins of the three atoms are affected homogeneously. The inhomogeneous evolution of 
the atoms causes decoherence, but this decoherence can be suppressed efficiently by applying a bias 
magnetic field of modest strength perpendicular to the plane of the atoms. The resulting lifetime of 
the RFF qubit can be many days, making RFF qubits of this kind promising candidates for quantum 
information storage units. Specifically, we examine the situation of three 6 Li atoms trapped in a 
CC>2-laser-generated optical lattice and find that, with conservatively estimated parameters, a stored 
qubit maintains a fidelity of 0.9999 for two hours. 

PACS numbers: 03.67.Pp, 03.65.Yz, 03.67.Lx 
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I. INTRODUCTION 

The processing of quantum information — be it for 
quantum communication, for quantum key distribution, 
or for quantum computation — requires the storage, ma- 
nipulation, and retrieval of qubits that are carried by 
physical systems. This "hardware" can consist of se- 
lected well-controllable degrees of freedom of single pho- 
tons, trapped ions, N-V centers in solids, or quantum 
dots, to name a few. All of them have advantages and 
disadvantages that make them well-fit for some purposes 
and unsuitable for others [0 . Practical applications that 
go beyond proof-of-principle experiments rely on qubits 
that are sufficiently robust for the task at hand. 

For example, fault-tolerant quantum computation re- 
quires a gate fidelity that is very close to unity. Lack 
of control over the system, however, always gives rise to 
decoherence. The typical decoherence time for quantum 
dots, ions in a trap, or diamond N-V centers is in the 
order of microseconds or milliseconds, and a decoherence 
time of a few seconds is in reach of the current technolo- 
gies [l|. This is still not quite sufficient for carrying out 
some complicated gate operations, nor for storage pur- 
poses 0. 

We explore here a scheme to overcome the decoher- 
ence problem with reference-frame-free (RFF) qubits 
constructed from three spin-1/2 neutral atoms. These 
RFF qubits have a remarkably long lifetime — a NMR 
proof-of-principle experiment that employs three spin- 
1/2 nuclei is on record [3| — and the alignment of ref- 
erence frames between observers, or the drift of frame 
between storage and read-out, is not an issue. 

The said construction of RFF qubits from trios of spin- 
1 /2 atoms was studied previously [H [j| . These atoms are 
individually highly sensitive to magnetic stray fields, but 
their symmetric RFF states are completely insensitive as 
long as the stray field affects all three atoms in the same 
way @. Decoherence of the RFF qubit could, therefore, 



result from spatial inhomogeneities of the stray field. We 
investigate the effect of such inhomogeneities and demon- 
strate that the RFF qubits are highly robust: For typical 
experimental parameters, the magnetic stray fields are of 
no concern. 

Rather, the lifetime of the RFF qubit is limited by the 
magnetic dipole-dipole interaction among the constituent 
atoms in conjunction with intrinsic imperfections of the 
experimental set-up. Our analysis, which uses conserva- 
tively estimated parameters and reasonable assumptions 
about experimental imperfections, shows that a stored 
qubit maintains a fidelity of 0.9999 or 0.999 for two or 
seven hours, respectively. 

The outline of the paper is as follows. We first review 
the construction of the RFF qubit in Sec. [TT1 followed 
by a qualitative discussion of possible causes of decoher- 
ence. In Sec. IIII1 we analyze the robustness of the RFF 
qubit against decoherence under random stray magnetic 
fields and conclude that the stray fields are innocuous. 
We then consider, in Sec. IIV1 the inhomogeneous mag- 
netic dipole fields of the partner atoms and find that, 
in view of unavoidable experimental imperfections, they 
are the dominating effect that limits the period for which 
quantum information can be stored. Section [V] briefly 
argues that the RFF system formed by three spin-1/2 
constituents is more robust than the system formed by 
four spin-1/2 constituents. An experimental realization 
could use ultracold atoms in an optical potential of a 
suitable geometry 0413 i we mention one possibility of 
creating a potential of this kind in Sec. lVIl We close with 
a Summary and Discussion. 



II. CONSTRUCTION OF THE RFF QUBIT 

For a system of three spin-1/2 particles, the total spin 
j of the system can be either 1/2 or 3/2. In the j = 3/2 
subspace, there are four states with different m values; in 
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the j = 1/2 subspace, there are two different m values, 
m = ±1/2, with two states each, which we distinguish 
by the quantum number A: \j — 1/2, m — ±1/2, A) with 
A = or A = 1. 

We construct these four orthogonal basis kets in the 
j = 1/2 sector by a variant of the procedure described in 
Ref . Q . When denoting the Pauli vector operator for the 
fcth atom by <X£, the total spin vector operator is given 
by 



1 3 

^E 



(1) 



k=l 



in units of %. The lowering operator J_ = J x — iJ y has 
the partner operators [Tl| 



Ok- 



(2) 



where a^- — {tJkx ~ iok y )/2 is the lowering operator for 
the fcth atom, and q — e l27r / 3 is the basic cubic root of 
unity. We write |±,A) for \j = 1/2, m = ±1/2, A) for 
brevity and choose the "+" kets in accordance with 
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and an application of J_ gives the corresponding 
kets, 

I- a) = j-\+,\) = j-Qx\m), 

with the outcomes 

l+,o) - (\in)q+\m)q 2 + \m))/V3, 
i+,i> = (im>9 2 + int>g + im»/V3 

and 

l-,o) = -{\m)q+\m)q 2 + \m))/V3, 
\-,i) = -{\m)q 2 + \m)q+\m))/V3. 
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The arrows symbolize "spin up" and "spin down" in 
the z-direction, so that |tft) = \j = 3/2, m = 1/2) in 
Eqs. ([3]) and ©. The |±, A) kets have the usual proper- 
ties as eigenstates of J z , namely 



J X |±,A) = |T,A)/2, 
J y \±,X) = |t,A)(±z/2), 
J Z \±,X) = |±,A)(±l/2), 



(7) 



as one verifies immediately 

Since the eigenvalues of J 2 distinguish the j = 1/2 and 
j = 3/2 sectors, we can express the respective projectors 
in terms of J 2 , 
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j=3/2 



(8) 



consistent with 1 
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J=l/2 



Pj=3/2 and J 
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Pj— 3/2- Alternatively, the projector onto the subspace 
with j = 1/2 is given by 

^■=1/2 = E E \^ x )(^m 

cr=± A=0,1 

= 7(3 - <J\ ■ er 2 - rr 2 ■ cr 3 - cr 3 ■ eri) , (9) 
6 

where the latter expression is available either as a conse- 
quence of Eqs. © and © or of Eqs. © and (JTJ) . 

The j = 1/2 subspace and the j = 3/2 subspace are 
both four-dimensional Hilbert spaces and, therefore, they 
can be regarded as tensor product spaces of two qubits, 
respectively. This is of no consequence for the j — 3/2 
sector, but it permits to write the j — 1/2 sector as com- 
posed of a rotationally invariant signal qubit — the RFF 
qubit — and an idler qubit [5|: the kets \<j, A) of Eqs. ©- 
© are labeled by the idler quantum number a = ± and 
the signal quantum number A = 0, 1. 

With the idler qubit in a maximally mixed state, the 
RFF state is identified by 

P RFF -E E W,X) 1 - Pxy (a,\>\= 1 -l 2 ®p^, 

cr=± A,A'=0,1 

. (10) 

where this tensor-product structure applies within the 
subspace with j = 1/2. We exhibit the statistical opera- 
tor of the signal qubit alone, 



-RFF 



\,\'=0,1 



|A)PAA'<A'| 



(11) 



by tracing over the idler qubit. The hermitian Pauli op- 
erators Ei, E2, E3 for the RFF qubit, 

i(Ei+zE 2 ) = |M)M|=l a ® (|0)<1|), 



<7 = ± 



£ 3 = E E mh-i)V,ai 

<y=± A=0,1 

= 1 2 ®(|0)(0|-|1)<1|), (12) 
are explicitly given by 



Ei + iT, 2 = -(<Ti • <r 2 + q 2 <J 2 ■ <J3 + q<?3 ■ o"i) ■ 
1 

<Ti • (<T 2 X <T 3 ) . 



(13) 



These are clearly rotationally invariant and possess the 
algebraic properties of Pauli spin operators in the j = 1/2 
subspace, such as (E1) 2 = Pj = \/2 and E1E2 = 1E3. 

Once the information is encoded in a RFF qubit (jlll) . 
with the idler in the maximally mixed state as in Eq. (|10p 
or in some other state, the information will be perfectly 
preserved as long as all three spin- 1/2 atoms precess in 
unison. In the non-ideal circumstances of a real exper- 
imental situation, however, the interaction with the en- 
vironment and the interactions among the physical car- 
riers of the qubit could cause decoherence, because the 
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atoms may be subject to torques of different strengths. 
Inevitably, there will be sources of noise over which the 
experimenter lacks control. It is our first objective to 
demonstrate that the information stored in the RFF 
qubit is preserved for a long time if magnetic stray fields 
with typical properties affect the carrier atoms. 



for any four vectors v%, v 2 , ^3, and V4 that pick out the 
components of G, whereby g 2 is the variance of the gaus- 
sian distribution of the diagonal entries of G, and T is 
the decay constant for the temporal correlation. We note 
that the diagonal elements and the off-diagonal elements 
of the gradient matrix do not have the same variance: 



III. MAGNETIC STRAY FIELDS 
A. Noise model 

The part of the Hamiltonian that describes the effect 
of the noisy magnetic field on the trio of atoms is given 
by 



fe=i 



(14) 



where /i B is the Bohr magneton (if necessary multiplied 
by a gyromagnetic ratio) and b k (t) is the randomly fluc- 
tuating magnetic stray field that acts on the fcth atom. 
The b k (t)s vanish on average, 



b k (t) = 0, 



(15) 



where the overline indicates the stochastic average. Since 
the atoms are close to each other, the fluctuations in 
the magnetic fields at the positions of the atoms are not 
independent but correlated. The dominant part of the 
noisy magnetic field is the same for all three atoms, and 
only a small part of the noise affects the atoms differently 
as a consequence of the nonzero gradient of the magnetic 
field. Upon denoting the gradient dyadic of the magnetic 
stray field by G(t), we have 



b k {t) - bt(t) = G(t) ■ (r k - n) 



(16) 



where r k is the position vector for the fcth atom, and 
G(t) is assumed to be independent of position within the 
small volume of relevance. This gradient component is 
the inhomogeneous noise that gives rise to decoherence of 
the RFF qubit, while the homogeneous noise is innocu- 
ous. 

In the noise model considered, every component of the 
homogeneous stray field and every component of the gra- 
dient dyadic has a random gaussian distribution with a 
vanishing mean. Owing to the Maxwell's equations, the 
gradient dyadic has to be symmetric and traceless: 



vi ■ G(t) ■ v 2 = v 2 - G(t) ■ vi , 
e a -G(t)-e o = 0. 



E 



(17) 



a—x,y,z 



It follows that the two-time correlation function of the 
field gradient has the form 



V! ■ G(t) -v 2 v 3 - G(i') • v A 

(3t>i • V 3 V 2 ■ V 4 + 3vi -V4V2-V3- 2vi ■ V 2 V 3 • v 4 ) 



xVg-rli'-t] 
4 y 



(18) 



e • G(i) • e e • G(t') • e = - e • G(t) ■ e' e • G(t') • e' 

3 



7 2 e -r|t'-t| 



(19) 



where e and e' are two orthogonal unit vectors. The 
autocorrelation function for the diagonal components of 
G is 4/3 times that of the off-diagonal components, while 
they have the same correlation time 1/r. 

As a consequence of Eq. (fT8|) . the autocorrelation of 
the field difference of Eq. (fl6| is 



Vx ■ (b k {t) - h(t)) (b k (t') - bl(t')) ■ V 2 



4^ 



3vx ■ v 2 r k 2 + vx ■ r H r ki ■ v 2 \ (20) 



with r k i — r k — r;. The correlation between the stray 
fields at the sites of the fcth atom and the Ith atom is 
then given by 



vi ■ b k (t) bi(t') -v 2 = e 



1 



b v\ ■ v 2 



-g (3vi-v 2 r ki + vi ■ r kl r kl ■ v 2 ) , (21) 

o 

where b 2 ^> g 2 r k \ 2 is the strength of the same-site corre- 
lation, as we recognize by a look at the k = I version, 



V!-b k (t) b k {t')-v 2 =b 



2„-r|t'-ti 



v 1 -v 2 . 



(22) 



Consistency with Eq. (|20|) is established by using Eq. ([21 
four times on the right-hand side of 

v 1 ■ (b k (t) - bi{t)) (b k (f) - b^f)) -v 2 



= vt ■ b k (t) b k {t') ■ v 2 + V! ■ bi(t) bi{t<) ■ v 2 



v x ■ b k {t) bi{t') • v 2 - v x ■ bi(t) b k {t') ■ v 2 . (23) 



Equations (jl5|) and (l2Tj) define the noise model that 
we use in Sec. IIII Cl below to derive the master equation 
by which we then study the effect of the inhomogeneous 
magnetic stray fields on the RFF qubit in Sec. IIII D[ and 
on a single-atom qubit in Sec. IIII El The model is char- 
acterized by the three parameters T, b 2 , and g 2 , which 
would have to be determined from experimental data 
when applying the model to an actual laboratory situ- 
ation. Other noise models are conceivable, in particular 
if one wants to describe a specific noise source of known 
characteristics. The noise model of Eqs. (fT5|) and ([21]) is 
generic, however, and quite suitable for the purpose at 
hand. 

In an experimental realization of the three-atom RFF 
qubit, there will be nearby Helmholtz coils for the precise 
control of the magnetic field at the location of the atoms. 
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-3/2 m=-V2 m=l/2 m=3/2 



228.2MHz 



/=3/2 



-1/2 m=l/2 



7=1/2 



2w 



FIG. 2: Ground-state hyperfine levels of the neutral 6 Li atom. 
The / = 3/2 quartet is separated from the f — 1/2 doublet 
by a transition frequency of 228.2 MHz. Three 6 Li atoms 
confined to their / = 1/2 ground states serve as the spin- 1/2 
particles from which the RFF qubit is constructed. 



FIG. 1: Three b Li atoms are trapped at the corners of an equi- 
lateral triangle. The probability clouds indicate the center- 
of-mass distributions whose spread w is about one-twelfth of 
a, the distance between the atoms. 



Typically, these coils are about 50 cm away and carry 
currents of about 1 A that are stabilized to 100 ppm or 
better [jjj]. Now, a current of 0.1m A at a distance of 
0.5m gives rise to a magnetic field of 4 x 10~ n T and 
a field gradient of 8 x 10~ n T/m. Not assuming any 
fortunate cancellation of the contributions from different 
coils, the values b = 5 x 10~ 10 T and g = 10~ 9 T/m are 
conservative estimates for these noise parameters. The 
temporal properties of the fluctuating currents tend to be 
dominated by the ubiquitous 50 Hz noise that the wires 
pick up, while high-frequency noise can be filtered out 
very efficiently, so that a correlation time of 1/r = 20 ms 
is a reasonable estimate [l2|. We will use these numbers 
throughout the paper. 



B. Lithium-6 

To be specific, but also mindful of possible experiments 
with two-dimensional confinement 13|, we consider the 
situation of Fig. [1] Three 6 Li atoms at the corners of 
an equilateral triangle, perhaps the sites of neighbor- 
ing minima of an optical potential such as the one dis- 
cussed below in Sec. |V] Each 6 Li atom is in the hy- 
perfine ground state with f — 1/2, which is energetically 
below the / = 3/2 hyperfine state by 2nh x 228.2 MHz; 
see Fig. H 

We denote the electronic spin operator of the feth atom 
by Sfc, with s k 2 = 3?i 2 /4, so that the energy of the atom 
trio in an external homogeneous magnetic bias field 6o is 
given by 



k=l 



(24) 



where we take the value of 2 for the gyromagnetic factor 
of the electron. With each atom confined to its / = 1/2 
ground state, this becomes 



H 



bias 



E 

k=l 



MB 



F k ■ b n 



(25) 



where —2/3 is the gyromagnetic ratio and Fk = (h/2)tJk 
is the total atomic angular momentum in the present 
context; the magnetic moment of the spin-1 nucleus is 
ignored. As indicated, here we identify the Pauli opera- 
tors of Eqs. ©, ©, (USD, or (HD. 

We choose the bias field in the z-direction, perpen- 
dicular to the xy-plane in which the atoms are located, 
bo = —Boe z , and express its strength in terms of the cir- 
cular frequency loq- ftuio — 2/i B _Bo/3; then 



1 3 

-ffbias = -fa^O ^ Okz = fu*)(>J z 
fc=l 



(26) 



The coupling of the f — 1/2 and the / = 3/2 multi- 
plets by the bias field is ignored, which is permissi- 
ble if the field is weak on the scale set by the energy 
difference, that is: w <C 2-7T x 228.2 MHz. For exam- 
ple, this condition is met for the modest field strength 
of B = 2mG = 2 X 10~ 7 T, when w = 2tt x 2 kHz 
is a thousandth of a percent of the transition fre- 
quency, and transition probabilities are of the order of 
(10- 5 ) 2 = 10" 10 . 

We need the bias field to fight the "internal magnetic 
pollution" that originates in the magnetic dipole-dipole 
interaction between the spin- 1/2 atoms. In terms of the 
electronic spin operators, this interaction energy is [lij 



2pb 

h 



(k,i) 



si - 3s fc • e k i e kl ■ sA 



(27) 

where a = \ri-i\ is the common distance between the 
atoms at the corners of the equilateral triangle, e k i — 
r k i/a is the unit vector that points from the fcth to the 
jth atom, and the summation is over the three pairs. As 
in the transition from Eq. to Eq. (|2"5j) . the restriction 
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Ml 



(x2) 



-■^frujo \ Ml 



,v""j = 3/2 



-§fiw + Ml 



and to transition frequency f2/(27r), which is 6mHz, if 
we continue to use the numbers found above. There is 
a clear separation of time scales, then, and noise with a 
correlation time 1 /Y in the order of 20 ms — which we 
regard as a typical number, see above — would not be 
able to induce transitions between the states separated 
by several Hloq while it will mix the states that are sepa- 
rated by Ml only or not at all. 

There are, of course, stray fields in the radio frequency 
range but their sources (the radio stations) are far away 
so that the gradient parameter g is extremely small, and 
noise of this kind is of no concern. By contrast, noise 
originating in nearby sources — current carrying wires 
in the vicinity of the laboratory, say — is relevant. 



FIG. 3: Level scheme for the effective three-atom Hamiltonian 
of Eq. (|32[) . The separations are not drawn to scale: tkoo is 
many orders of magnitude larger than hQ. The two j = 1/2 
levels are degenerate doublets; this degeneracy is exploited 
for the encoding of the robust signal qubit. 



1 „ h 

Sk -> -r-Pfc — 

3 6 



to the / = 1 /2 ground state amounts to the replacement 

(28) 

which turns Eq. (l27l) into 

#dd = \jiSl ^2 ( a k ■ vi - 3(T k ■ e M e H ■ <r{) (29) 



(k,i) 



with 



Ml 



Mo Mb 
4tt 3a 3 



(30) 



For a distance of a = 883 nm (see Sec. lVII below). we have 
fl = 2tt x 6mHz, smaller than ujq by a factor of 3 x 10 5 , 
so that the transitions induced by H d d are completely 
suppressed in the presence of a 2mG bias field. There- 
fore, only the part of Hm that commutes with i/bias of 
Eq. (|26|) is relevant, and we arrive at 



H, 



1 



dd 



Ml (3 j\ 



(31) 



as the effective Hamiltonian for the magnetic dipolc- 
dipole interaction among the atoms. Note that this H^d 
vanishes in the j = 1/2 sector where the signal and idler 
qubits reside. 

The combined effective Hamiltonian 



H dd — hw J 



hn(3J* - j 2 ) 



(32) 



has the non-degenerate eigenvalues zL^Hloq + Ml and 
±^Mjq — Ml in the j = 3/2 sector, and the two-fold 
eigenvalues ±ifiw m the j = 1/2 sector; see Fig. [31 The 
energy differences correspond to transition frequencies of 
about 1, 2, 3 x loq/(2it), which are in the few-kHz range, 



C. Master equation 

In view of this separation of time scales — very fast 
wo-oscillations and very slow 51-oscillations on the scale 
set by the correlation time 1 /Y of the random stray field 
— we can use master equation techniques to account for 
the net effect of the noise. For the purpose of deriving 
the Lindblad operators of the master equation, we put 
i/dd of Eq. (|32[) aside and use an interaction picture in 
which the fast wo-oscillations of -Hbias are transformed 
away. The statistical operator in this interaction picture, 
denoted by g(t), then obeys the von Neumann equation 
of motion 



d_ 

df 



W) - T [Q(t), flnoisoW] 



(33) 



with 



xt . (f\ _ iH hias t/hrr . (f\ p -iH hias t/H 
iJ noise \ L J — c L± noise 



iojQtJ ' _ 



'(-f )!>(*) 



• cr k e 



-iujQtJ z 



(34) 



fc=i 



where the replacement fj, B — > —fj, B /3 accounts for the gy- 
romagnetic ratio that we first met in the transition from 
Eq. (HID to Eq. (US]). 

The unitary evolution operator U(T) links g(T) to the 
initial statistical operator g(0), 



g(T) = U(T)g(0)U(Ty . 
We solve the Lippmann-Schwinger equation 

U{T) = l- l -J^ dtH noisc (t)U(t) 
to second order in H nn ^„, 



(35) 



(36) 



U(T) ~ 1 - iMT) - ^4>i{T? ~ iMT) , (37) 
where the hermitian phases </>i(T') and <p2 (T) are given 
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by 



and the stochastic averaging of Sec. IIIlAI turns this into 



2ih 2 J 



dt H noise (^) 



dt / dt' [H noisc (t),H Doisc (t')] . (38) 



To second order in H no [ SC , then, we have 



g(T) = g(0)+i[e(0),MT) + <h(T)] 



p{T) = p(0)+i\p{0), MT)] 

+ Umt)M0),Mt)]] 



(40) 



where p(t) — g(t) and the initial statistical operator is 
not affected by the averaging or the transition to the in- 
teraction picture: g(0) = g(Q) = p(0) = p(0). Note that 



(39) 0i (T) = follows from Eq. flS 



We take a closer look at the "sandwich term" in the double commutator, 

3 r T r T 



MT)p(0)MT) = (g 



■J pi pi 

J2 / dt / dt'a k {t)- b k (t)p(0)b t (t') -atf). 
k,i=i Jo J ° 



(41) 



Here, T is much longer than the correlation time of the noise, TT ^> 1, so that there are very many cycles of the uiq- 
oscillation in a short ^-interval. Therefore, the rapidly oscillating terms in cr k {i) do not contribute to the ^-integration 
and the replacement 



Vkz^z 



is permissible; and likewise for <x;(i'). This "rotating- wave approximation" takes us to 



MT)p(0)MT) 



g0 <r kz p{0)ai z dt dt' e z -b k (t) b t {t' 

u i i JO JO 



)-e z 



k,l=l 
3 



= (i) E^(°k 



; ( ff a) 2 (l-^) 



fc,i=l 



T /-T 

/ 

o Jo 



dt / dt'e- 1 " 1 *-*' 1 



(42) 



(43) 



after using Eq. ([2~Tj) for t>i = t>2 = e 2 , r^; • e z = 0, and with the Lindblad operator £ given by 



(1 — 8ki)a . For TT ^ 1, the remaining double 



integral equals 2T/T, and we arrive at 



Mt)p(o)Mt) 



" ^2<7kzP(0)<Jkz + ^-jJzP(0)Jz 



4r 



k=l 



with the time constants 



3r 



p B ga 



and 



3M 2 



86 2 -3( 5 a) 2 



(44) 



(45) 



(46) 



Since ga -C b, we have t' -C t, and the numbers of 
Sec. UTTBl that is: 1/T = 20 ms and ga = 9 x 10~ 16 T, 
give t — 2 x 10 10 s — an amazingly long time. 

The replacement of Eq. (|42[) gives a vanishing commu- 
tator in the double integral for 02 (T) in Eq. (|3"5|) . so that 



4>2(T) = in Eq. (|40[) . In summary, then, we have 

p(T) ~ p(0) + T£p(0) (47) 



1 3 1 



fe=i 



and the master equation in the interaction picture is sim- 
ply 



^p(t) = Cp(t) 



(49) 



Upon getting out of the interaction picture, and re- 
introducing the slow ri-oscillations of _£/dd, this gives us 
the master equation 

§£p(t) = jj\P®, ^bias + H dd ] + Cp(t) (50) 

for the evolution of the coarse-grain, stochastically aver- 
aged, statistical operator p(t). 

We note in passing that this master equation could al- 
ternatively be derived with standard textbook methods, 
such as those that proceed from the Redfield equation, 
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here: 



_ ?(T) = _ / dt [ Hn 



,XT),[g(t),H noise (t)}] , (51) 



and invoke the Born-Markov approximation and the 
rotating-wave approximation to arrive at Eq. (1491). For 
details of this procedure, see section 3.2 in Ref. [HI], for 
example. 

In view of the diagonal form of the Lindblad operator 
in Eq. (|4"5|) . the expectation value (A) t = tr {Ap(t)} of 
an observable A obeys the differential equation 



dt 



(A) t 



in 



H dd ]) t + (CA) t 



(52) 



For observables that commute with J z , which is the case 
for all operators related to the signal qubit, only the J 2 
part of .f/dd and the r-term of C matter, and then the 
simpler equation 



d 

di' 



(A) t = f([A,J 2 ]) t + l-j^{a kz Aa kz ) t - ^(A) 



fe=i 



applies. In particular we have 

(f(J z )) t = {f(J z )) 
for all functions of J z . The product 
Ho hT 



Or = 



-5 



47T a 3 (ga) 2 



(53) 
(54) 

(55) 



states the relative size of the time constants in Eq. (|53j) . 
It depends rather strongly on the distance a between the 
atoms; we have fir ~ 10 9 for the values used earlier 
(a = 883 nm, 1/r = 20 ms, ga = 9 x KT 16 T). 



D. Time dependence of RFF-qubit variables 



By making use of Eq. (|53|) , we can now calculate 
(Pj = i/ 2 ) t , the probability of finding the three-atom sys- 
tem in the j = 1/2 sector at time t, and the time- 
dependent expectation values of the RFF-qubit Pauli op- 
erators of Eqs. (fl"3")) and (fT"2"|) . The outcome is 



,-th 



(Ei + *E2> t = e-i*/ r <£i + »E2> 0! 
<E 3 ) t = e - t/r <£ 3 ) , 



(56) 



if the system is initially in the j = 1/2 sector, where 
the idler and signal qubits reside. The expressions for 
(Pj=i/ 2 ). and (£3)4 are exact solutions of the respec- 



tive versions of Eq. (1331 . and the result for (Ei + «E 2 ) t 
is an approximation that neglects terms of relative size 



0.75 




0.25 



FIG. 4: Comparison of the data from a numerical simulation 
with the analytical results of Eqs. (|56[) . Curve "a" displays 
(Pj = i/2) t , curves "b" show (£i) t for (£i) = 1, and curve 
"c" is for (Sa) t with (£3)0 = 1. The crosses are from a 
simulation of the dynamics, averaged over 1000 runs. The 
solid- line curves represent the analytical results of Eqs. (|56|) . 
The dotted "b" curve shows what one would get for (£i) t if 
Q vanished rather than being large on the scale set by r; we 
observe that the inter-atomic dipole-dipole interaction accel- 
erates the decay of (Si).. 



Supporting evidence is provided by numerical simula- 
tions These are done by generating a RFF state 
at initial time t = and letting the state evolve with 
stochastic random noisy magnetic fields on the three 
atoms. All components of the noise are generated ran- 
domly at each time step from a gaussian distribution, 
and the Maxwell's equations are strictly imposed on the 
magnetic fields. Figure |4] shows both the analytical and 
the numerical results of the evolution of the RFF qubit 
for an arbitrary initial RFF state. We note that there is 
very good agreement between the results of the simula- 
tion and the analytical solution of the master equation. 

Quantum information stored in the RFF qubit is de- 
graded substantially only after a good fraction of r has 
elapsed. But since t — 2 x 10 10 s is more than 600 years, 
we conclude that the effect of the inhomogeneous mag- 
netic stray fields is of absolutely no concern. Put differ- 
ently, the experimenter need not take special measures 
to suppress the stray fields. 



E. Decoherence of a single-atom qubit 

If — rather than making good use of the three-atom 
RFF signal qubit — one encoded quantum information 
into the / = 1/2 ground state of a single 6 Li atom, the 
effect of the random magnetic stray field would be de- 
scribed by the single-atom master equation 



d_ 

di 



Pit) 



IUJQ 



1 



W),o z ] + —[a z Mt)^A] (57) 
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with 

Ti = — = -At) T = T \2^b) <<T - (58) 

The resulting time-dependent expectation values are 

{a x +ia y ) t = e^e-^i^+iay)^ 

Wt = Wo. ( 59 ) 

so that the quantum information can be stored for a frac- 
tion of time t\ . 

The bias field stabilizes the z component: It sepa- 
rates the spin-up and spin-down states in energy by Hujq 
and so prevents transitions between them — this is, of 
course, the essence of the rotating-wave approximation 
of Eq. ((32). Therefore, one could encode a classical bit 
in a single spin- 1/2 atom and protect it from the stray 
magnetic field. |17j | 

Without the bias field, the master equation 

§tP® = -^["-Mt),<r]] (60) 
applies. Its solution 

p(t) = ±(l + e- 2t ^{*) -cr) (61) 

shows that the state decays toward the completely mixed 
state with a life time of ri/2. It follows that, in addition 
to preserving the z component, the bias field also slows 
down the decay of the x and y components by a factor 
of two. 

For the example used in Sec. IIII Al — a fluctuating 
0.1 mA current in a wire at a distance of 50 cm — we have 
ga/b ~ 10~ 6 and obtain t\ ~ 2 x 10~ 12 t. Even for the 
very large value of r found above, r ~ 2 x 10 10 s, the life- 
time of the single- atom qubit is quite short: t± ~ 40 ms. 
Clearly, the well-protected RFF qubit of the three-atom 
system has an advantage over the unprotected single- 
atom qubit: The stray fields, which are of no concern for 
the RFF qubit, have a devastating effect on the single- 
atom qubit. 

More relevant than the lifetime t\ is the duration of the 
initial period of high fidelity. We recall that the fidelity 
of two single-qubit states, specified by their respective 
Pauli vectors, is given by 



F(pi,p 2 ) = ^(l + s 1 -s 2 ) + ^l-s\^l-sl. 

(62) 

For the fidelity F(t) = F(p(t),p(Q)) between the initial 
qubit state and the state at later time t, we have the 
lower bound 




so that a fidelity of, say, 0.999 is only guaranteed for a 
fraction of a millisecond. By contrast, the RFF qubit 
would have a fidelity of 0.9999 or better for several 
months if nothing mattered except for the magnetic stray 
field. 

We note that the ratio of ri and r is solely determined 
by the comparison of the distance between the atoms and 
the distance of the atoms from the source of the noise, for 
which we have been using 883 nm and 50 cm, respectively, 
implying ga/b = 900 x 10~ 9 /50x 10 -2 ~ 2xl0~ 6 . There- 
fore, the conclusion that n/r ~(2x 10~ 6 ) 2 ~ 10~ 12 
holds irrespective of the actual physical process that gen- 
erates the stray magnetic field as long as the noise source 
is half a meter away. 



IV. DIPOLE-DIPOLE INTERACTION 

Two assumptions of ideal geometry enter the deriva- 
tion of the effective Hamiltonian for the dipole-dipole in- 
teraction in Eq. (|3~Tj) : That the atoms are located at the 
corners of a perfect equilateral triangle; and that the bias 
field is exactly perpendicular to the plane of the atoms. 
Let us now consider the consequences of imperfections on 
both counts. 



A. Center-of-mass probability distribution 

As illustrated by the probability clouds in Fig. [TJ the 
atoms do not have definite positions but rather prob- 
ability distributions for their centers of mass, given by 
the ground-state wave functions of the trapping poten- 
tials. We assume that, for the purpose at hand, the re- 
spective trapping potentials are reasonably well approxi- 
mated by isotropic harmonic oscillator potentials, so that 
each atom has a gaussian probability distribution, 

p(r) = (V2^wy V2 e-^ r2 / w2 , (64) 

where r = is the position of the trap center and w 
is the width of the gaussian. The oscillator frequency 
wtra P /(27r) of the trap is related to w and the mass M of 
the atom by 

Wtrap = 2ih ' (65) 

which is obtained by fitting the potential around the bot- 
tom of the trap to a harmonic-oscillator potential. For 
the potential of Sec. I VII below, Wtrap — 2tt x 0.3 MHz, 
so that the width of the gaussian wave function is w = 
75 nm. Accordingly, here, earlier in Fig. [TJ and in what 
follows, we take the width w to be about one-twelfth of 
the distance a between the atoms. 

When comparing the restoring force of the oscillator 
potential, Muj t rap 2 r, with the dipole forces exerted by 
the partner atoms, Ml/ a, we find that the balance of 
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forces would shift the equilibrium position by an amount 
of the order of 



29. 



Muj tlap 2 a 



Wtrap a 



4 X 10~ 9 W, 



(66) 



which is a completely negligible effect. We also note that, 
depending on the joint spin state of the three atoms, 
the shift is in different directions, and the center-of-mass 
degrees of freedom get entangled with the spin degrees 
of freedom but, since the shift is such a tiny fraction of 
the position spread w, this entanglement is so weak that 
it can be safely ignored. As a consequence, the center- 
of-mass motion is decoupled from the dynamics of the 
spins, and probability distributions as in Eq. ([64]) apply 
to the atoms at all times. 

The total statistical operator for the three-atom sys- 
tem is then the product g(t)g CM of the spin factor g(t) 
of Sec. IIII Cl and a static center-of-mass factor g CM . The 
von Neumann equation for g(t) is obtained by tracing 
over the center-of-mass variables, 

d % 

g- t Q(t) = ~ [e(t), tr CM {g CM H tot }] , (67) 



where H tot = H cu 



H nois C 1S the t0tal 



dd bias 

Hamiltonian. Of its four terms, the dipole-dipole inter- 
action energy H dd and the noise part H noise involve both 
spin variables and center-of-mass variables. In view of 
the lesson learned in Sec. IIII1 however, there is no need 
to deal with H noise in detail. 

We consider the center-of-mass average of the contri- 
bution from atoms 1 and 2 to i?dd, 



tr, . 



^CM^dd J 



Mo / Mb , _ _ 

4tt V 3 ' 



x / (dn)(dr a ) 



Pi(r 1 )p 2 (r 2 ) 



ri + r 2 



(68) 

where a is the vector from the trap center for atom 1 to 
the trap center for atom 2. Allowing for different widths 
of the two gaussians, the integration yields 



(dn)(dr 2 ) 



Pi(ri)p 2 (r 2 ) _ }_J 



a-ri + r 2 \ a \ ^/2w\ + 2i 



(69) 

with the standard error function erf( ). Its asymptotic 
form 



erf(z) = 1 - 



• for z > 1 (70) 



tells us that the right hand side of Eq. differs from 
1/a by a term of relative size 10 -16 for W\ ~ w 2 — a/12. 
It follows that 



, / rr(12)\ MO /Mb 



-v a v a - 

a 



&2 



(71) 



in the present context, which is exactly the (k,l) — (1,2) 
term in Eq. (f2l?|). and the center-of-mass probability dis- 
tribution is of no further concern. 



B. Non-ideal geometry 

But we need to account for the unavoidable imper- 
fections of any experimental realization: The triangle 
formed by the trap centers for the tree atoms is not ex- 
actly equilateral, and the plane of the actual triangle is 
not exactly perpendicular to the z-axis defined by the 
bias field. 

First, with au denoting the distance between the fcth 
and the /th trap center, we define the average distance a 
by means of 



E 

(k,i) 



(- 



(72) 



and measure the deviation of au from a by the small 
parameter au, 



( — J =l-3a H , Va fc( =0. (73) 

(k,l) 



This average a value is used in Eq. (|30|) to determine the 
dipole-dipole coupling strength Ml 1 and we have 

Hdd = iftO J^(l - 3a k i)cr k ■ (l - 3e M e fe; ) • er ; (74) 

(M) 



instead of Eq. (|29|) . where 1 denotes the unit dyadic. If 
the atoms are indeed trapped in the minima of an optical 
potential, \a k i \ — 10 -2 is achievable without resorting to 
extreme measures [l~8j . 

Second, nonzero z-components of the unit vectors e kl 
require 

1 - 3e a e M -> i(3e z e z - 1) [l - 3(e 2 • e kl f] (75) 



for the step from Eq. ([29]) to Eq. ([31]) . Misalignments 
that exceed 1° can be avoided with standard experimen- 
tal techniques ;18|, so that \e z ■ e k i\ = 10~ 2 is a conser- 
vative estimate. 

The combined effect of both imperfections is a modifi- 
cation of -ffdd, such that 

Hdd = (3J Z 2 - j 2 ) + nriK (76) 

o 

with 

K = e k i~(o- k ■ a i - 3(7 fcz g; z ) , (77) 
(M) 

rather than the e k i = version of Eq. (]3~T]) . The relative 
size of the imperfections is measured by 

(ki = o-ki + (e z ■ e ki ) 2 - 3a k i{e z ■ e u ) 2 , (78) 
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wherein, for the values of \au\ and \e z ■ e k i\ above, the 
three terms are of the order 10~ 2 , 10~ 4 , and 10~ 6 , re- 
spectively, and the a k i contribution dominates. 

We note in passing that the imperfection parameters 
aki and the dot products e z -e^i that appear on the right- 
hand side of Eq. ([78]) are not independent of each other. 
Rather, the triangle condition a\2e\2 +023623 +031631 = 
imposes the restriction 



^(l-3a fci )" 1/3 e 2 -e kl =0, 



(79) 



(fe,0 



where the summation over the pairs is cyclic, that is: 
(fc,/) = (l,2), (2,3), and (3,1). 

The operator K vanishes in the j = 1/2 sector, 



^ = 1/2^ = 1/2=0, 



(80) 



and there are no first-order contributions from the K- 
term to the evolution of the RFF qubit. It follows that, 
during the initial period of high fidelity, the geometrical 
imperfections contribute in second-order of the small e/y 
parameters. 



C. Time dependence of RFF-qubit variables 

For a quantitative analysis, we employ the master 
equation that results when Eq. (|50l) is modified in ac- 
cordance with the observations made here and above, 



d_ 

df 



P {t) = two !>(*), J,] + iV[p(t), Jz 2 - \J 2 + K] 



n 



[j z ,[p(t),Jz]}, 



(81) 



where we put 1/r — > 0, 2r' — > T\ in the Lindblad op- 
erator of Eq. (|48|) . While the double-commutator term 
leads to the fast decay of single- atom spin coherence, as 
we saw above in Sec. MI El it is of no consequence for the 
RFF qubit because all RFF observables as well as their 
commutators with J 2 commute with J z . It follows fur- 
ther that the expectation value (A) t of a RFF variable 
A, with the three-atom system initially prepared in the 
j = 1/2 sector, is given by 



(A) t = (UesityAUesit) 



(82) 



with 



U eB (t) = P J=1/2 exp(-int(J z 2 - § J 2 + K)) P j=1/2 
1 + /(*) 



-P 

2 

WW 



3=1/2 

T,i cos ip — E 2 sin 1^9) , (83) 



where 



/(*) 



fixe 



-iQ,2 1 



fil + fi 2 



9. 2 



} = ^(V(l-e) 2 +8| K | 2 ±(l-e) 



(84) 



for 



£ = £12 + £23 + £31 > 0, 



£12 + q £23 + 0£3i = e v n 



(85) 



For |a fcJ | ~ \e z ■ e kl \ ~ 10~ 2 in Eq. $7%$, we have 
first £ ~ \k\ 2 ~ 10~ 4 and then Oi ~ ft and Q 2 = 
2\K\ 2 n 2 /fl 1 ~ 10~ 4 f). Accordingly, /(*) is the sum 
of a long-period oscillation with a large amplitude and 
a short-period oscillation with a small amplitude; see 
Fig. El 

The resulting time-dependent expectation values of the 
RFF observables are 



( P i=l/2) t = 



1 



\f(t)\ 2 l-\f(t)\ 2 /^ 

2 2 X^ 1 /o 



l/(*)l a / s <v0 



!/t = /(^E^-iE 3/Q 



with 



Ej cos<y9 — E 2 sin^ , 



E 2 = S 2 cos (/? + E x sin (/? . 



(86) 



(87) 



Equations (|86j) contain all information about the state of 
the signal qubit in the course of time. We use them to 
evaluate the purity of the signal-qubit state at time t and 
its fidelity with the initial signal-qubit state. 

The purity [1 + s(t) 2 ]/2 of the signal-qubit state is 
quantified by the squared length of its Pauli vector, 

s(t) 2 = (^) 2 t+(^) 2 t +(^) 2 t 



= 1 - 



(^=l/2>* 

l/WI 2 



[i-«(o)1, 



so that an initially pure signal-qubit state, s(0) = 1, re- 
mains pure. When the initial state is mixed, s(0) < 1, 
both s(t) > s(0) and s(t) < s(0) are possible, depending 
on the relative size of \f(t)\ and (Pj=i/ 2 ) t - Specifically, 
we have 



s(t)^s(0) if (1 + |/W|)('l + (E^)J 



^2 (89) 



and |/(t)| < 1. 

The fidelity F(t) of the signal-qubit state at the later 
time t with the initial state is given by 



i-Hgn 



F(ty - , 

4(P J =i /2 ) t 
\f(t)\- Refjt) 
2<^=i/ 2 ) t 
It is bounded from below by 



>- S (0) 2 



F(t)>Jl- 



IWffll 2 

;i + i/wd 



2 ' 



(90) 



(91) 
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F(t) 



0.99995 



0.9999 




D. Compensating for triangle distortions 



15 30 

fiif/(27r) 




In Sees. IIVBI and IIV CI we regarded the imperfection 
parameters a k i and e z ■ e k i as resulting from the lack 
of perfect control over the apparatus, and their values 
would not be known with high precision. Suppose, how- 
ever, that the experimenter has diagnosed the set-up and 
knows the actual shape of the triangle quite well while 
having very precise control over the direction of the mag- 
netic bias field. She can then attempt to adjust the bias 
field such that the three efc/s of Eq. (JTSJ) are equal, with 
the consequence that k = in Eq. (|5o| and fit) = 1 in 
Eq. (|84p . We do not discuss this matter in further detail 
and are content with mentioning that, for small values of 
the afeiS, a bias-field direction e z with 



\ 



IP 

(k,l) 



Oikl 



(93) 



achieves this, where the approximation neglects terms of 
second and higher order in the a k is. With this compen- 
sation for the imperfections in the shape of the triangle 
by a judicious tilt of the bias field, the ratio ft 2 /fti can 
be reduced by much, with a corresponding lengthening 
of the initial period of high fidelity. 



o 1 1 * 1 1 

2500 5000 7500 10000 

«ii/(27r) 

FIG. 5: Fidelity of the RFF qubit. For fi 2 /fii = 10~ 4 , 
the plots show F(t) of Eq. I|90p and its lower bound of 
Eq. (J9T) for t < 45 x 27r/fti (top plot), for t < 150 x 2?r/rii 
(inset in the top plot), and for t < 2n/Q.2 (bottom plot). 
Curve 'a' is the lower bound on F(t); the other three curves 
are for (£>[ v) } = 0.4 and s(0) = 1 (curve 'b'), s(0) = 0.8 
(curve 'c'), and s(0) — 0.6 (curve 'd'). One can clearly see 
the small-amplitude short-period oscillations and the large- 
amplitude long-period oscillation. For the parameter values 
used throughout the paper, the respective time ranges are 2, 
7, and 450 hours. 

where the equal sign holds for f(t) ^ 1 if, for example, 
s(0) = 1 and (Z[ v) ) = 2/(1 + \f(t)\) - 1. For il 2 < Ox 
and t <§; 27r/f22, this bound is 

n 2 t n 2 . n 2 t 

F(t) > cos — - — sin — + • • ■ (92) 

where the ellipsis stands for terms of order {Vt 2 /Vti) 2 . 
This two-term approximation serves all practical pur- 
poses for f^/^i — 10~ 4 . The fidelity is assuredly very 
high during the early period dominated by the small- 
amplitude oscillations with frequency fii/(27r) ~ £1/ (27r): 
We have F = 0.9999 or better for 45 periods of the fast 
Six oscillations, and F = 0.999 or better for 140 periods, 
when Vl 2 = 10 — 4 Sli . These matters are illustrated in 
Fig. El 



V. RFF QUBIT FROM FOUR SPIN-1/2 ATOMS 

An alternative construction of a RFF qubit uses four 
spin-1/2 atoms and their two-dimensional subspace with 
j = 0. A pure state of the RFF qubit is then realized by 
a pure state of the four-atom system, and decay results 
from leakage to the sectors with j — 1 and j = 2, which 
are nine-dimensional and five-dimensional, respectively. 
By contrast, a pure state of the three-spin- 1/2-atom 
RFF qubit corresponds to a mixed state of the three- 
atom system with leakage into the space of the idler qubit 
and into the j = 3/2 sector. Clearly, the two construc- 
tions of the RFF qubit are substantially different. 

For comparison with the two-dimensional equilateral 
triangular configuration of three atoms, let us consider 
four atoms located at the corners of a square; see 
Fig. MX)- The system is stabilized with a bias magnetic 
field perpendicular to the plane of atoms to reduce the 
decoherence due to the internal pollution from the dipole- 
dipole interactions. The distance between the two diago- 
nal pairs of atoms is V2 times larger than the distance be- 
tween the four pairs at the sides. Thus, the dipole-dipole 
interaction is unbalanced between all pairs and it cannot 
be made rotationally invariant by the bias magnetic field. 
Explicitly, the effective dipole-dipole interaction is here 
given by 

H d d = ^Kl^iJl - J 2 - 8 ^ [3(fTi 2 a 3z +a 2z a 4z ) 
- (erj • cr 3 + cr 2 • cr 4 )] ) (94) 
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(i) 3 (ii) 

FIG. 6: RFF qubit constructed from four spin-1/2 atoms. 
(i)Two-dimensional square configuration. The dipole-dipole 
interaction is unavoidably unbalanced here because the dis- 
tance for the two diagonal pairs is larger than the distance for 
the four edge pairs, (ii) Three-dimensional pyramidal config- 
uration. Here, if the height h is chosen such that b/a = 0.661, 
the effective dipole-dipole interaction has equal strength for 
all six pairs of atoms. 




nt/(2vr) 

FIG. 7: Lower bounds for the fidelity F(t) (solid curve) and 
the expectation value {Pj=o{t)) (dashed curve) for the RFF 
qubit constructed from four spin-1/2 atoms. 

and f(t) is exactly of the form in Eq. (|84p with fii and 
2 replaced by 



with Ml as in Eq. (|30|) . For the same reasons as in the 
three-atom case, the stray magnetic field is of no concern, 
and the system evolves unitarily 

p(t) = e -i(H a +H dA )t/h poel (H a + H dA )t/h (gg) 

The projector onto the j = subspace of the four-atom 
RFF qubit is given by 



Pj=o — 7;{Si2S34 + + S14S23) ! (96) 



where the SjkS are the singlet state between jth and fcth 
constituents, 



(97) 



for j, k — 1,2,3,4 and j 7^ k. 

The effective dipole-dipole Hamiltonian has the struc- 
ture of Eq. ([To]) , with the operator K now given by 

K = j(<Tl • <T3 + <T 2 • CT 4 - 3(T lz a3 z - 3<72zCT4z) , (98) 

where c = (4 — \/2) /6 accounts for the relative reduction 
in the strength of the dipole-dipole interaction for the 
two diagonal atom pairs. Moreover, Eq. (1521 continues 
to apply for the expectation value of a RFF operator A, 
with the four-atom system initially prepared in the j = 
sector. Thereby, the effective unitary evolution operator 
is now 

U eS (t) = P J=0 e- lH ^ t/h P 3 = 

= i± /^o + + V3E 2 ) , (99) 

where S^, £ 2 , and £ 3 are the RFF Pauli operators for 
the four-atom system as defined in Eqs. (25) of Ref. [5j, 



fix 
Q 2 



n 



= 'j (v/(2- C ) 2 



8c 2 ± (2 - c) 



( 1.780 
\ 0.2m 

(100) 

5, very different from 



which have a ratio of r2i/f2 2 : 
the ratio of 10 4 in Eqs. ([84]) and 

The expectation values of the RFF operators in the 
j = sector can be obtained analytically, in particular 
for the expectation value (Pj=o) t and the RFF qubit fi- 
delity F(t), for which the obvious analogs of (-P/=i/2/ t in 
Eq. (HH) and F(t) in Eq. ([90]) apply. With the four-atom 
version of /(£), the lower bound on F(t) of Eq. (|9T|) is 
valid, and we also have (Pj = o) t > \f{t)\ 2 . Both lower 
bounds are shown in Fig. [7] for the high-fidelity period 
of < fit < 2tt x 0.112. We have a fidelity of 0.9999 
or better for t < 0.062 x 2tt/0 and 0.999 or better for 
t < 0.11 x 2tt/CI. 

Note, in particular, the substantial probability of los- 
ing the four- atom RFF qubit: After the lapse of t = 
0.08 7 x 27r/f2, there is a chance of more than 10% that 
the four-atom system has left the j = sector. This 
is a consequence of the rather small f2i/T2 2 ratio. By 
contrast, for the three-atom qubit with Oi/f2 2 = 10 4 , 
the persistence probability \Pj=i/2/ t is never less than 
0.9996. 

The dipole-dipole coupling strength f2 is proportional 
to 1/a 3 , where a is the length of the sides of the square. 
If we use laser beams with the same wavelengths as 
used for the three-atom system in Sec. I VII to construct 
the potential for four atoms in a square geometry, the 
inter-atomic distance is a = 663 nm and Eq. (|30j) gives 
n = 2tt x 16 mHz. It follows that we can guarantee 
F > 0.9999 for about four seconds and F > 0.999 for 
about seven seconds. This shows that even with a per- 
fect square geometry, the RFF state constructed from 
four spin-1/2 atoms decays about 2000 times faster than 
the qubit constructed from three spin-1/2 atoms in an 
imperfect equilateral triangle configuration. 

The imbalance in the dipole-dipole interaction strength 
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FIG. 8: Six coplanar laser beams consist of two sets of three 
coherent beams; the angle between beams within each set is 
2tt/3. The respective wave vectors have lengths |fcj.| = \k%\ ~ 
\ks\ = 2ir/\ and \k^\ = |fes| = |fcs| = 2n/X'. Different lattice 
structures can be created by alternating the phases of the 
laser beams. For the lattice of our design, we keep the set of 
beams with wavelength A to be in phase and the phases for 
the other three beams with wavelength A' are 2n/3, 0, and 
-2tt/3. 



between the pairs of atoms could be removed by using the 
three-dimensional pyramidal configuration of Fig. [BJii). 
To fight against the internal magnetic pollution we ap- 
ply a bias magnetic field perpendicular to the plane of 
three of the atoms which are arranged in an equilateral 
triangle, with the fourth atom above the center of the tri- 
angle. By adjusting the height of the pyramid such that 
the ratio between the atomic distances is b/a — 0.661, a 
rotationally invariant effective dipole-dipole potential is 
obtained. The lifetime of the RFF qubit for this four- 
atom geometry is comparable to the lifetime of the RFF 
qubit constructed from three atoms in the equilateral- 
triangle configuration examined in the previous sections. 
But the practical realization of the peculiar pyramidal ar- 
rangement, a distorted tetrahedron with reduced height, 
is much more challenging than the equilateral triangle. 

Clearly, there is no advantage in using the RFF qubits 
made from four atoms over the RFF qubits made from 
three atoms. Rather, the simpler three-atom system is 
preferable. 



VI. STRUCTURE OF THE OPTICAL LATTICE 

One kind of optical lattices that could be used is a mod- 
ification of the Kagome lattice, where we have an equi- 
lateral triangular lattice of which every site is formed by 
three spin- 1/2 atoms arranged in a small equilateral tri- 
angle. A possible physical construction of such a lattice 
is to use two sets of three coplanar coherent laser beams 
and the angle between the beams within each coherent set 
is 27r/3. By arranging them in the configuration shown 
in Fig. [8] and adjusting the phases, an optical trapping 
potential with the contour plot shown in Fig. [9] can be 
produced. For the potential presented in Fig. [9l we chose 
to keep the phases of the set of beams with the longer 
wavelength A to be the same and the phases for the other 
three beams with the shorter wavelength A' are 27r/3, 0, 
and — 27r/3, but this choice of phases is not unique. 

Under the joint consideration of being able to address 




FIG. 9: The top figure shows the contour plot of the overall 
optical potential, and the three global minima at each lattice 
site of the big triangular lattice indicate where the trio of 
atoms can be trapped. The bottom left figure gives a more 
detailed view of the site where one trio is trapped. The dashed 
curve in the bottom right figure is a plot of the potential 
along the vertical line cutting through the global minima a 
and b; and the solid curve is a plot of the potential along 
the horizontal line cutting through the other global minima 
c. The two cuts intersect at the saddle point s. 



each RFF qubit individually, of the necessity of high laser 
intensity to trap the 6 Li atoms, and of the requirement 
of a low probability that the atoms scatter photons from 
the trapping lasers, we propose to use the CO2 laser, 
from which the desired wavelengths can be generated by 
frequency doubling (or tripling) in nonlinear media. CO2 
lasers with wavelength 10.6/^m are widely used and a 
high beam intensity is routinely achieved, and optical 
trapping of 6 Li is reported [l9j]. We can thus have A = 
10.6 /im and A' = A/8 = 1.33 /im for the wavelengths 
in Fig. [5J For the produced lattice the atomic distance 
between atoms within the trio forming one RFF qubit is 
a — 883nm and the distance between the RFF qubits is 
7.1/xm. The potential strength associated with the low- 
frequency laser is four times that of the high-frequency 
laser for the example of the potential given in this section. 

The recoil temperature of the D\ (or D2) line of 6 Li 
is T rec = 3.5 fiK. To have an estimation of the trapping 
frequency, the function of the potential at the minima is 
fitted with a spherically symmetric harmonic potential. 
With the ratio between the laser intensity and the re- 
duced saturation intensity Io/I s = 10 s , the frequency 
of the harmonic potential is about w tr ap — 1.8 MHz. 
The energy separation of the first excited state and 
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the ground state of the harmonic potential is given by 
Huj trap ~ 1.8 x 10~ 28 J (=13 fiK), which is about five 
times the recoil energy. The polarizability of 6 Li is 
24.3 x 10 _24 cm 3 , which yields a scattering rate below 
10 _3 s~ 1 at this intensity, corresponding to a scattering 
time of more than 1000 s for one photon per atom. Con- 
sequently, the recoil heating is negligible. Upon assigning 
a gaussian profile to the center-of-mass wave function of 
the atoms, the spread of the wave function is estimated 
to be w = 75 nm for the trapping frequency w trap above. 

The depth of the optical trap is roughly 1.1 x 10~ 27 J, 
corresponding to a temperature of about 80 /iK and more 
than twenty times the recoil energy. The required inten- 
sity of the 10.6 /an CO2 laser is about Jo = 0.2W//zm 2 
and for the 1.33 /im laser it is about I' = 0.03 W//im 2 . 

The numbers given above for the trap are for the sam- 
ple potential presented here, which serves the purpose 
of demonstrating that such a lattice can be had. The 
properties of the trapping potential depend on the laser 
intensities Iq and 7 . Other methods for making an array 
of atomic trios are conceivable. This is a hardware issue 
and all details are determined by the experimental set-up 
at hand. 



VII. SUMMARY AND DISCUSSION 
A. Summary 

We studied the effect of stochastic magnetic stray fields 
on the RFF qubit made from three spin- 1/2 atoms and 
found that the RFF qubit decoheres very slowly although 
the spin states of the individual atoms decay very quickly. 
The only coupling of the spin- 1/2 atoms to the environ- 
ment is through their magnetic dipole moments, so that 
magnetic stray fields give rise to uncontrolled changes of 
the quantum state of the atoms. The long lifetime of 
the RFF qubit results from its insensitivity to the over- 
all magnetic field and its fluctuations because they affect 
all three atoms equally and, therefore, do not affect the 
RFF qubit at all. Decoherence of the RFF qubit origi- 
nates in spatial variations of the magnetic field, but they 
are subject to the constraints imposed by the Maxwell's 
equations. For parameter values that are typical for ex- 
perimental situations, we find that the RFF qubit can 
maintain a very high fidelity for months — if the mag- 
netic stray field is the only source of decoherence. 

We then analyzed the effect of the dipole-dipole in- 
teractions among the three atoms and imperfections in 
the geometry of the trapped atoms. We found that the 
inter-atomic interactions bring more decoherence to the 
RFF qubit than the fluctuating magnetic stray fields, al- 
though the dipole-dipole interaction itself is a unitary 
process. Nevertheless, the RFF qubit states were shown 
to be very robust within the parameter regime and un- 
der our assumptions. As an example, we showed that the 
period of high fidelity, say F = 0.9999, with respect to 
the initial state is roughly two hours. 



B. Assumptions 

Let us review the assumptions that we used to derive 
the result and discuss their validity. 

We have assumed that the atoms are trapped in the 
deep optical lattice at very low temperature so that their 
center-of-mass motions are negligible. One example for 
such a desired optical lattice, created by standard laser 
techniques, was presented in Sec. ED We estimated, in 
Sec. I IV Al the effect of the center-of-mass motions and 
found that it is much smaller (i.e., 10~ 16 ) than the effect 
of the dipole-dipole interactions for the optical lattice 
considered. 

The most challenging element seems to be to maintain 
the stable lasers for the optical lattice in order to ob- 
serve the long-time evolution of the RFF qubit. In a real 
experiment the collisions with rest-gas atoms are also in- 
evitable and they may very well limit the lifetime of the 
RFF qubit in practice. We note that drifts of the lasers 
in time would not spoil the long lifetime as long as all 
lasers are locked in phase. This is because the time scale 
of these parameter changes is much slower and hence all 
atoms follow the optical lattice adiabatically. 

The noise model employed in this study is divided into 
two types. One is a fluctuating magnetic stray field aris- 
ing from unavoidable imperfections in the surrounding 
apparatuses such as the Helmholtz coils, electric wires, 
and so on. The Helmholtz coils used to generate the ho- 
mogenous bias magnetic field are identified as the major 
source for the noise of this kind. Other possible fluc- 
tuating magnetic fields are much smaller than this and 
less inhomogeneous as the respective sources are farther 
away. We then linearized these fluctuating fields around 
the homogeneous bias field to analyze the decoherence 
for the RFF qubit in Sees. HEl and HTCl 

The other type of noise is due the magnetic dipole- 
dipole interaction among the three atoms. We have 
shown that these inter-atomic interactions are a major 
source of decoherence for the RFF qubit when imperfec- 
tions of the experimental set-up are taken into account. 
In Sees . lIVBl and ITVQ we accounted for deviations from 
the ideal equilateral triangle configuration of the three 
atoms as well as a misalignment of the magnetic bias 
field and found that such insufficiencies still allow for a 
very long lifetime of the RFF qubit. In fact, we observed 
that imperfections in the geometry of the three atoms 
could be compensated for by adjusting the direction of 
the bias magnetic field, provided that the experimenter 
has sufficient control over the relevant parameters. 

We should not forget to mention that the rotating- wave 
approximation was used, for example, when analyzing 
the effect of fluctuating magnetic stray fields. This ap- 
proximation is valid within the energy scale of our set-up, 
where a probability for a non-resonant transition is very 
small. The numerical study without the rotating-wave 
approximation also supports the validity of our master- 
equation analysis. 

All these results are, of course, derived with the as- 
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sumption that the parameters used in this study are well 
controlled with a certain precision. We however took 
rather conservative numbers for these parameters so that 
we can estimate a realistic lifetime of the RFF qubit. 
This also leaves some room for improving the lifetime of 
the RFF qubit in further studies. 



advantages and disadvantages depending upon the phys- 
ical system one chooses. All these are largely unexplored 
territories that need further surveying. 



D. Outlook 



C. Alternatives 

In the scheme presented here, the RFF state is made 
from three spin- 1/2 atoms. But this is not the only 
way to construct RFF qubits. We can equally well use 
three identical atoms of any non-zero ground-state spin j, 
where the RFF signal qubit lives in the subspace of total 
angular momentum 3j — 1, which has two states for each 
m value. The idler space is then (6j — l)-dimensional. 
Analogously, RFF qutrits can be constructed in the sec- 
tor of total angular momentum 4j — 1 from four iden- 
tical atoms with ground state spin j. Such alternative 
constructions offer considerable flexibility in choosing the 
isotope for a practical implementation. 

As we mentioned in Sec. [Vj other geometries and four 
atoms are not better than the case of three atoms in 
the equilateral triangle, but they are still good for the 
purpose of storing quantum information. For example, 
three or four atoms on a line, with the bias held in the 
right direction, could be an easier choice for a trapped- 
ion experiment. Likewise, physical systems other than 
cold atoms should in principle give a long lifetime if in- 
formation is encoded in the RFF subsystem. There are 



Having thus shown how to construct robust units for 
storing quantum information, the problems of how to en- 
code and read out the information, and how to process 
it in such a set-up, need to be addressed. We will report 
progress on this front in due course. 
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